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A 

a,D 

2c 

ce2(n,-q) 

E 


=  (cosh2$  -  cos2n) 

=  serai-major  and  semi-minor  axes  of  the 
ellipse  with  £  =  £0 

=  interfocal  distance  of  ellipse 

=  Mathieu  function  of  the  first  kind,  of 
the  order  of  2,  defined  by  Equation  (A1.6) 

=  Young's  modulus 

=  eccertricity  of  ellipse 


Fek-  =  Fek,(£,-q)  = 


(Fek2)  ' 


(Fek2)" 


H,B,C,D,E ,F,M 
hl,h2 

In(z)'Kn(z) 


J 

k 

l 

M 

P 

(Po)  ' 


=  k-type  modified  Mathieu  function  of  the 
second  kind,  of  the  order  of  2,  which  is 
defined  by  Equation  (A1.4) 

=  J^[Fek2(5,-q)] 

a2 

=  — 2[Fek2  (E,  ,-q)  ] 

=  F/ (c/2) 2 

=  superposition  constants 

=  scale  factors 

=  modified  Bessel  functions  of  the  first  and 
second  kinds,  respectively,  where  n  de¬ 
notes  or.der 

=  the  Jacobian  of  transformation 


1/2 


k  £ 

2  i 


couple-stress  characteristic  length 


=  M- 


(P2)  ' 

=  uniaxial  uniform  load 


=_ce2(0,q)  co2  {J,q)/;.^2) 


|*  X  * **** 


se2  =  se2(n,-q) 


(se2)  • 


(se2)  " 


y  ,y 
x,py 


VyB 


V% 


a  ,o  .0  ,o 
xx  xy  yx  yy 


=  1^-1^]  =  Here  q  is  the  real  posi¬ 

tive  -humber 


=  semi-major  and  semi -minor  axes  of  any 
ellipse 

=  polar  coordinates 

=  Mathieu  function  of  tha  first  kind,  of  the 
order  of  2,  which  is  defined  by  Equation 
(A1.5) 

=  ^tse2(n,-q> 

a2 

=  — rtse0  ( n ,-q ) 
an2  2 

=  Airy  stress  function 

=  Cartesian  coordinates 

=  orthogonal  curvilinear  coordinates 

=  defined  by  Equation  (3.35) 

=  defined  by  Equation  (AIII.17) 

=  the  inclination  of  the  curve  B  -  a  constant 
to  the  x-axis  in  orthogonal  curvilinear 
coordinates  (a,B) 

=  modulus  of  rigidity 

=  couple-stress  compc  *ents  in  Cartesian 
coordinates 

=  couple-stress  components  in  orthogonal 
curvilinear  coordinates 

=  couple-stress  components  in  elliptical 
coordinates 


=  Poisson's  ratio 


=  ke 


=  ke 


=  elliptical  coordinates 

=  stress  comDonents  in  Cartesian  coordinates 


IV 


WV  -»>J 


■a  iMimu 


aa'  a8'  3a'  88 


=  stress  components  in  orthogonal  curvi¬ 
linear  coordinates 


a«'aen'V'°nn 


*  stress  components  in  elliptical  coordi 
nates 

=  coup.le-stress  function 
=  3ocal  rigid  rotation 
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1 . 0  INTRODUCTION 


1. 1  Objective  of  Research 

The  objective  of  this  research  is  to  obtain  solutions 
for  an  infinite  tension  p..ate  bounded  at  the  interior  by  an 
elliptical  hole.  The  nominal  tension  in  the  plate  is  uniform 
along  the  majo^.  axis  (see  Fig.  1)  .  The  couple-stress  effect 
is  considered.  Two  limiting  cases  for  the  problem  are  for 
the  ellipticity  of  0  and  1  for  the  interior  boundaiy.  That 
is,  for  the  first  case  the  interior  boundary  is  a  circle,  for 
the  second  case,  a  crack  along  the  major  axis.  Of  course,  the 
problem  of  a  crack  is  the  major  interest  of  the  research.  The 
present  investigation  includes  a  study  of  a  class  of  solutions 
which  satisfies  the  static  equilibrating  traction  on  the 


Figure  1. 
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1.2.  Review  of  Pertinent  Literature 


The  origin  of  couple-stress  theory  of  elasticity  is 
attributed  to  E.  and  F.  Cosserat,  with  modern  developments  by 
Truesdell,  Toupin,  Grioli,  and  Aero  and  Kuvshinskii.  A  dis¬ 
cussion  of  the  origin  and  development  of  the  couple-stress 
theory  is  given  by  Mindlin  and  Tiers ten  [1] . 

Recently,  specific  plane  problems  in  couple-stress  elas¬ 
ticity  have  been  studied  by  many  investigators.  The  following 

is  a  list  of  those  studies  which  are  related  to  the  present 

/ 

problem.  Mindlin  [2]  found  the  couple-  stress  effects  on  the 
stress  concentration  factors  for  a  circular  bole  m  various 
two-dimensional  fields  of  stress  when  the  diameter  of  the  hole 
is  comparable  in  size  to  the  couple-stress  constant  l. 

Weitsman  [3],  [41  generalized  the  solution  to  the  cases  of  cyl¬ 
indrical  inclusion  in  fields  of  cylindrical  symmetric,  end  om- 
axial  tension,  respectively.  Sternberg  [5],  [bj  found  the 
effect  of  couple-stresses  on  the  stress  concentration  around  a 
crack  assuming  that  the  stress  singularities  at  the  ctac..- 
tips  are  of  the  same  order  as  those  in  the  well-known  classi¬ 
cal  solutions.  An  earlier  report  by  Ju  and  Hsu  [7J  (AFOtT-o9- 
1908  TR)  contains  a  comprehensive  review  of  the  theory  and  the 
basic  equations  for  the  plane  problems.  All  pertinent  nota¬ 
tions  and  definitions,  therefore,  will  be  referred  to  [7]  to 
avoid  repetition. 


2.0  TWO-DIMENSIONAL  COUPLE-STRESS  THEORY 


The  approach  toward  the  solution  of  a  crack  problem  may 
be  either  by  means  of  the  degeneration  of  an  elliptical  hole 
[7]  or  by  the  use  of  a  half-plane  [5] .  The  choice  of  using  an 
elliptical  hole  has  the  advantages  of  (a)  no  a  priori  assump¬ 
tion,  (b)  ready  check  for  the  problem  of  a  circular  hole  [2] , 
(c)  ready  check  for  the  classical  solution  [8] .  In  order  that 
the  boundary  may  be  described  by  a  coordinate  line,  the  ellip¬ 
tical  coordinate  is  used.  Hence,  al]  equations  for  the 
couple-stress  theory  of  elasticity  are  expressed  in  such  co¬ 
ordinates  in  this  problem. 

2.1.  Rectangular  Cartesian  Coordinates 

According  to  Ju  and  Hsu  [7]*  or  Mindlin  [2],  the  field 
equations  of  two-dimensional  couple-stress  theory  in  a  state 
of  plane  strain  without  body  forces  or  body  couples  are  given 
as  follows: 

U  and  'F  must  satisfy 

V4U  =  0  (2.1) 

and 

V2 ( T  -  £2V2T)  =  0  (2.2) 

respectively,  where  l  is  the  couple-stress  characteristic 

2 

length  so  that  l  is  the  ratio  of  the  Cosserat  modulus  to  the 

5 

Numbers  m  brackets  designate  references  at  end  of  the 
report.  Equations  in  Section  2.1  refer  to  Equations  (1.40) 
through  (1.51)  m  [7], 
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modulus  of  rigidity  (y)  [7] .  U  and  V  must  also  be  related  to 
each  other  by  the  Cauchy-Riemann  equations 


JL(¥  -  £2V24')  =  -2(1  -  v)i2~(V2U) 

dX  oy 

-  l2V2V)  =  2(1  -  V)£2  ^r(V2U) 


(2.3) 


Stresses,  couple-stresses,  displacements,  and  rigid-body  rota¬ 
tion  are  expressed  in  terms  of  the  two  generating  stress  func¬ 
tions  U  and  ¥  as  follows: 


2  2 
d  U  _  By 

°xx  ,,  2  9x9y 


2  2 

°yy  3x2  ^3y 


?  9 

_  £U  _  33 

xy  “  9x9y  2 
3y 

(2.4) 

2  2 
9  U  9  Y 

°yx  9x9y  +  3x2 


94' 

9x 


My  '  9y 


2.2.  Orthogonal  Curvilinear  Coordinates 

To  generalize  the  results  [2],  [7]  which  are  derived  for 
rectangular  Cartesian  and  polar  coordinates,  consider  those 


5 

* 

f 

i 

f 


I 


]! 


for  the  general  two-dimensional  orthogonal  curvilinear  co¬ 
ordinates  (a, 8).  The  two  generating  stress  functions  U(a,0) 

and  4'(a,0)  must  satisfy  (2.1)*  and  (2.2),  respectively.  Here 

2  4 

the  Laplacian  operator,  V  ,  and  the  biharmonic  operators,  V  , 

are  expressed  in  terms  of  orthogonal  curvilinear  coordinates 

a  and  0. 


and 


ih  . 

h 

V2  =  h.h„ 

3 

i  i- '  +  i- 

2  8 

12 

8a 

h2  3a  j  36 

hl  3&, 

(2.5) 


V  = 


2  2 
V  V 


where  h^  and  h ^  are  scale  factors  in  o  and  0,  respectively. 
The  Cauchy-Riemann  equations  (2.3)  can  be  transformed  to  a 
and  0  coordinates  by  use  of  the  chain  rule  that 


8 

8a 

8  , 

38 

8 

33?  " 

8x 

8a 

8x 

38 

0  _ 

8a 

3  + 

30 

8 

8y 

8y 

8a 

3y 

38 

The  stresses,  couple-stresses,  displacements,  and  rigid-body 
rotation  are  expressed  m  terms  of  U  and  y  as  follows: 


Numbers  in  parentheses  are  references  to  equations  m 
the  text. 
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=  h 


act 


2  lil 

2  382 


^i.^lU+h  ^  i2  .  h  h  1?3L 

h2  3a  3a  2  3g  3g  12  3a38 


h  ii  h  ^  ii 

2  3g  3a  "  1  3a  38 


88 


2 

.2  32U  .  ^1  3U  ^2  ^1  3U  .  .  32V 

1  1  3a  3a  "  h1  38  38  1  2  3a38 


.  h  ^1  il  +  u  ^2  3f 
n2  38  3a  1  3a  38 


a8 


3hl  3U  ,  ^2  3U 

l“2  3a3 8  “2  38  3a  “  ftl  3a  36 


,  u  32U  _  H 

TTTn"  n 


2 

h2  ill  a-  'll  Ull  il  _  h  1^1  32 

n2  2  '  h2  3a  3a  ”  2  38  38 


(2.1 1 


8a 


h  h  ^  h  d-^L  i£  h  —L  i£ 

12  3a3 8  “  2  38  3a  1  3a  36 


2 

■  h2  iii  +  h  il  ^2  ^1  32 
nl  ,2  ni  3a  3a  '  h,  36  36 

O'*  X 


,  3  V 

U  =  h,  — - 

a  1  oa 


-  h  11 

u6  2  38 


2.3.  Elliptical  Coordinates 

The  elliptical  coordinates  £  and  n  (see  Fig.  2)  are  re¬ 
lated  to  the  rectangular  cartesian  coordinates  x  and  y  by 


6 


x  =  c  cosh£  cosn 


(()<_£;<<».  0  ^  n  <  ?  7T )  (2.8) 

y  =  c  sinh£  sinn 


Figure  2. 


Here  the  curves  £  =  constant  and  n  =  constant  form  an  orthog¬ 
onal  system  of  confocal  ellipses  and  hyperbolas,  with  the  com 
mon  foci  being  the  points  (±c,0.  Figure  3  shows  two  limiting 
cases  when  the  elliptjcity  e  is  1  and  0,  respectively. 

The  Jacobian  of  transformation  is 


At  5  =  0,  n  =  0  or  it,  J  =  0.  Hence  the  only  singular  points 

of  the  transformation  are  located  at  the  foci  (±c,0) .  The 

scale  factors  hr  and  h  can  be  calculated  by 

£  n  J 


ah5  _  ahn  _  ft  ”2 
TT"  35  ~  c 


sinh2£ 


ft  "2  .  , 

—  A  sin2n 

L> 


h2  -  2  a"1 

h  —  — A 

n  7 


/3f  **2 

h  =  h  =  —  A  Z 

5  n  c 


3h  3h  _ 

h£  7|”  =  h£>  a|“  =  2  A  sinh2- 

C 


(2.11) 


n  3n  "  5  3n 


— j  A  sin2ri 

c 


Equation  (2.5;  becomes,  by  use  of  (2.9)  and  (2.10), 


V2  =  4 


4  +  i 

9  5  cm 
L 


(2.12) 


4  2  2 

V*U  =  V*\TU 


The  Cauchy-Riemann  equations  (2.3)  can  be  transformed  to 
elliptical  coordinates  5  and  n  by  use  of  (2.6)  with  a  =  5, 
(i  =  0 


3  (^  _  5, 2 v 2 T )  =  -2(1  -  vU2  ^-(V2U) 

dt,  dn 


-  i2V2\<)  =  2(1  -  sj)l2  |^IV2U ) 


(2.13) 
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rum 


Similarly,  stresses,  couple -stresses,  displacements,  and 
rigid-body  rotation  are  expressed  in  terms  of  elliptical 
coordinates  by  use  of  (2.7)  and  (2.11). 


_  2  n-l  3  U  ^  2  .-2  .  .  3U  2  .-2  _  3U 

Op,  —  — j  A  — j  +  — 2  ^  sinh25  -jr-r  —  — ^  A  sin2n  — — 

^  c  3n  c  **  c  dn 


2  .-1  3  4*  2  .-2  .  0  3  4* 

— ar  A  t*tt —  A  sin2n 

c2  3^3n  c2  3? 


A"2  smh25 
c2  an 


(2.14) 


2  -1  3  0  2  -2  .  ..  dU  2  -2  „  dU 

i  =  — =•  A  — -  -  — ~  A  smh25  -r— -  — =■  A  sm2n  — 

nn  c2  c2  35  c2  3n 


+  i_  A“1 
c2  A  353n 


2  ,-2  0  3Y 

2  A  sin2n  3  £ 


A  2  Sinh25 
c2  an 


(2.15) 


5n 


2  ,-l  3  U  2  -2  0  3U  2  -2  .  3U 

7  A  353^  +  ^  A  Sin2n  35  +  ^  A  Sinh2^  cTn 


2  -i  ai 
T 

c 


3  2y  2  .-2  „  3Y  ^  2  .  -2  -  3f 

— ■“  — 7T  A  smh2£  tt*  +  — ??  A  s in 2 n  - — 

3n2  c2  ^  c2  on 

i 


2.1i 


n5 


2  A-1  32U  2  .-2  _  3U  2  A-2  u«,  dU 

— *r  A  —  +  — sp  A  Sin2n  rryr  +  — sr  A  Sinh25  — — 

c2  353n  c2  35  c2  s  dn 


+  2  -i  ri 

c*  35 


Ur  = 


(2) 


1 

'2  ay 

35 


A-2  sinh25  ~  \  A-2  sin2n 

c2  c2  9n 


(2.17) 


(2.18) 


u„  = 


-  1 

(2) 2  "2  3y 

c  3  n 


(2.19) 
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Components  of  stress  and  couple-stress  in  elliptical  coordi¬ 
nates  are  shown  in  Pig.  4. 


y 


2.4.  Boundary  Values  of  the  Problem 

The  couple-stress  solutions  are  obtained  for  the  uniform 
tension  plate  bounded  at  the  interior  by  an  elliptical  hole. 
The  nominal  tension  p  is  parallel  to  the  major  axis.  For  the 
present  general  class  of  problems,  the  boundary  conditions  at 
the  elliptical  hole  are 


o)^(n)dn  =  0 


n=n 


(2.20) 
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LIUW.UI 


1 


E=5CJ 

n=n 


u>2(n)dn  =  0 


f2ir 

i 


u2(n)dn  =  o 


n=n 


(2.21) 


(2.22) 


where*  a>1(n)  and  w2(n)  are  some  weiqhting  functions  defined  on 
the  interior  boundary. 

The  regularity  conditions  at  infinity  are 


o°°  =  p ,  a°°  =  a 00  =  a°°  =  u00  =  u00  =  0 

yy  c  xx  xy  yx  x  y 


(2.23) 


3.0  SOLUTION  OP  THE  PROBLEM 

3.1.  Selection  of  Stress  Functions  U(g,n)  and  ¥(g,ri) 

2 

As  U(£,n)  is  biharmonic  in  (2.1),  with  7  taking  the 
form  of  (2.12),  we  choose  five  solutions  of  U(£,n) 

=  e2^  +  cos2n 
-2£ 

U2  =  e  +  cos2n 

U3  =  e~2^  cos2n  (3.1) 

U4  =  « 

U5  =  e2^  cos2n 

1 

The  selection  of  these  fundamental  biharmonic  functions  is 
based  on  the  same  argument  as  that  used  by  Filon  and  Coker  [8], 
Notice  that  (2.2)  can  be  rewritten  as 

V2  ('F  -  £2V24')  =  0  (3.2) 

and 

(1  -  £2V2 )  V24'  =  0  (3.3) 

From  (3.2),  we  obtain  the  wave-function  solution  as  the 
products  of  two  types  of  Mathieu  function  (see  Appendix  I). 
Here  we  choose 

4*2  =  (Fe*2)  (se2)  (3.4) 
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Equation  (3.3)  implies  that  one  solution  of  S'  is  harmonic,  S'1« 

furthermore,  Cauchy - Ri amann  equations  (2.13)  conclude  that 

2  2 

this  harmonic  S'  and  2(1  -  v)£  V  u  must  be  conjugate  harmonic 
functions  [2]  [7],  As  such,  S^  is  selected  with  U  in  (3.1), 
as 


uf  _  sin2  n  /  o  c  \ 

1  ~  (cosh2£  ~  cos2n) 

Finally,  0  and  S'  are  formulated  by  a  superposition  of  stress 
functions  given  in  (3.1),  (3.4),  and  (3.5). 

U  =  HTl^  +  BU2  +  CU3  +  DU4  +  EU5  (3.6) 

S'  =  FS>1  +  MS<2  (3.7> 


3.2.  General  Expressions  of  Stresses  and  Couple-Stresses 

Define  o_P  ,  o  ,  o  ,  a  ,  \i~  ,  and  u  equivalent 
^i  nni  <’,‘i  n*i  ^i  ni 

to  IK  for  i  =  1,2, 3, 4,  and  5  and  S*^_3j  for  i=6  and  7-  Thus 

we  have 


0 =  HO  r  r  +  Bo  _  _  +  C  Orr  +  Do  _  r  +  E  0Pr  +  FoPP  +  MoPP 

«  u5  ZK6 


a„  =  He;  +  Bam„  +  Co„„  +  Da  +  Eo  +  Fa  +  Ma 
nn  nni  nn2  nn3  nn4  nn5  nn6  nn? 


o  P  =  Ha,  +  Bo  P  +  Ca  >  +  Da  _  +  Eo  ,  +  Fa„P  +  Ma  P 

n  ih2  r)€3  n?4  nC5  nS6  n?7 


aP  =  Ho.  +  Ba  _  +  Ccr  +  DaP  +  EoP  +  FaP  +  MaP 

Sn  Cn1  Cn2  5n3  £n4  5n5  Sng  $n? 
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For  U2#  we  have 


1  -2 

are  =  — A  [2  cos4n  -  8  cos2n  cosh2£ 
^ 2  c2 


+  4  +  2e“4^] 


o__  =  -ij-  A  2  [2  cos4n  -  8  cos2ne  2^  +  4  +  2e  4^] 


nn 


2  c 


1  -2 

or  ~  °  c  -  — y  A  [-4  sin2n  cosh2£l 
^n2  n^2  c 


v2u2  =  —  A-1[e“2?  -  cos2n] 


y-V2u2  -  — jjj-  A  2  sin2n  sinh2£ 
c 


-^-V2U2  =  -  A  2[1  -  cos2n  cosh2£] 


/ (3.10) 


For  U-j,  we  have 


af{.  =  Ar  A  ‘[003411  •  e  -  cos2n  (e  +  3)  +  3e  2^]' 
”3  c‘ 


2  ,-2 -  .  -2€  _ 


a  =  -tt  A  2[-cos4n  *  e  2^  -  3e  2^  +  cos2n(e  4’  +  3)] 
n3 


a_  =  a  -  =  A  2[sin4n  •  e  2^  -  sin2n  (e  +  3)] 

Cn3  n^3  c2 


2”  -  3  V2U„  =  A-  V2'j,  =  0 


(3.11) 


yc3  ~  yn3  V  u3  3n  ’  "3  n  ’  J3 
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.1  ■  «■'  '".W 


For  U^,  we  have 


acr  =  -V  A-2  [2  sinh2£l 
c 


a  =  -4  A-2 [-2  sinh2£l 

nrijl  2 

4  C 


a-  =  a  r 

5n4  nS 


■  A  2 [sin2n3 
4  c 


Mr  -  »„  =  V204  =  ±  V2U4  -  |f  V2t)4  =  0 

4  4 


(3.12) 


For  U5,  we  have 


orr  =  -4-  A  2 [2  cos4n  •  e2^  -  cos2n(2e4^  +  6)  +  6e2^]  \ 


«5  c2 


a  =  -4- A-2  [-2  cos4n  *  e2^  +  cos2n(2e4^  -t  6)  -  6e2^] 


nn5  c2 


ar  =  a 


7  V 


=  -4- A  2  [-2  sin4ri  •  e2^  +  sin2n(2e4^  +  6)] 


/ (3.13) 


£nc  nCc  2 


ur  =  v  =  v2u,  =  2-  v2u5  =  ^  v2u5  •«  0 

n5  5  3ri  5  dt,  5 


For  4'1,  we  have 


=  4r  A~4  sinh2£[l  +  sin22n  -  cos2n  cosh2£] 


££ - 1 

”6  c 


a  =  A-4  sinh2Ul  +  sin22n  -  cos2n  cosh2?] 

m6  7 
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TJ/JP  ■'•MW  MI  *.j'  *' ■  ■_ 


Of.  =  a  _  =  A  4  sin2n [sinh22£  -  1 

^n6  n^6  c 


+  cos2n  cosh2£] 


i  =  -  A  [1  -  cos2n  cosh2£] 
n6  c 


2/2  -5/2  r  .  _  .  . 

ur  =  -  — —  A  [sm2n  sinh2?] 

^6  c 


(3.14) 


=  0 


4’1  -  S,2V24'1  =  A_1  sin2n 


~(V 
3  n  1  1 


2  2 

fi.  V  ¥1) 


-2 

-2A  tl  -  cos2n  cosh2^] 


— n* 
hi 


2  2 

2,  V  4^) 


=  -2A  2(sin2ri  sinh2£] 


For  ^2 '  we  have 


aCS  =  •  -^"2  [-A(Fek2)  ’  (se2)  1  +  (Fek2)  ’  (se2)  sin2n  ' 


+  (Fek2) (sinh2C) (se2) ’ ] 


nn< 


— j[-A(Fek2) ' (se2) '  +  (Fek2) ' (se2)  sin2n 


c  A 


+  (Fek2) (sinh2?) (se2) ' ] 
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0 E,r\  =  ~2~2  t_A<Fek2^se2^  "  '  (sinh2^)  (Fek2)  '  (se2 ) 
7  c  A  i 


+  U‘'ek2Mse2)'  sin2n] 


>(3.15) 


r  £  =  ~ 2~2 ^A^Fek2) " (ge2^  ~  (sinh2£) (Fek2 ) ' (se2 > 

**7  cA* 


+  (Fek2)(se2)1  sin2n] 


=  I'D 


1/2 


h  c 


IaI  (Fek2)'(se2) 


1/2 


.1/2 

\  =  $\z\  (Fek2)(se2)> 


*2  “  “  *VV 


^<*2  *  =  0 


} 


Substituting  the  expressions  in  (3.9)  to  (3.15)  into 
(2.14)  to  (2.19),  the  results  are  as  follows: 


2  4  f 

=  2  2^H^cos‘*'1  ”  ^  cos2r]  cosh2£  +  2  re  ■') 

c  A 


+  B[cos4v  -  4  cos2n  *  e2^  +  2+  e^D 
+  C[cos4n  •  e  2^  -  cos2n(e  ^  +  3)  +  3e  2D 


+  D[sinh2£] 
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tg**n4J 


•VJ«>s*r  -'se^^s'  >W'Wj>  ^»~vr*"r  -  _  -  -,._  „  -  .  yj  ,  -  ^  .  *-  -  ‘  ■  '*' 

+  E[cos4n  •  e2^  -  cos2n(e4^  +  3)  +  3e2^] 

-  F[-^-  sinh2£(l  +  sin22n  -  cos2o  cosh2£)J 

.  A1 

•■  M[-A(Fek2) ' (se2) '  +  (Fek2 ) ' (se2 ) sin2n 
+  (Fek2) sinh2C (se2) ' ] }  (3.16) 

a  =  _ L_{H[cOS4ri  -  4  cos2rie2^  +  2  +  e4^] 
nn  C2A2 

-2  ?  -4£, 

+  B[cos4n  -  4  cos2ne  ^  +  2  +  e  ] 

-  C[cos4n  •  e  2^  +  3e  2^  -  cos2n  (e  4^  +  3] 

-  D [sinh2EJ 

-  E  [cos4n  •  e2C  +  3e2^  -  cos2n(e4ii  +  3)] 

+  F[-%  sinh2Ul  +  sin22n  -  cos2n  cosh2C)l 
A2 

-  M[-A(Fek2) ' (se2) '  +  (Fek2) ’  (se2) sin2n 

+  (Fek2) sinh2? (se2)  '  ] )  (3.17) 

a , n  =  2  2{ 2  U-  -  B)  [sin2n  cosh2£] 

c  A 

+  Clsin4n  •  e-2^  -  sin2n  (e  4^  +  3) ] 

+  D[sin2r|] 
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fit* wsr**K. 


+  M[  -A{Fekn)  (se7)  "  -  sinh2£  (Fek0)  '  (se7) 


+  (Fek2) (se2) '  sin2n]} 


’n £  ~~  %  +  ~ 2-  M[A(Fek2P(se2)  +  A(Fek2)  (se2)  "] 
c  IK 


y  -  — 1.2  {_2£  i&j  n2ii  sinh2fj  +  M[A2  (Fek2)  '  (se2)  ] ) 

^  cA  / 


(3.18) 


(3.19) 


(3.20) 


t/o  ^  o 

un  ~  • ---5-^2 ^ -2F  [1  -  cos2n  cosh2£j  +  M [A  (Fek2 )  (se2)' ]  )  (3.21) 
cA 


I  The  seven  unknown  cocf  ;i  chonts  ’(,  S,  C,  D<  K ,  7,  and  M  in 

(3.15)  to  (3.21)  will  be  determined  by  using  boundary  condi¬ 
tions*  at  £  =  £q,  Cauchy-  -Riemann  equations,  and  regularity 
!  conditions  at  i  y  ^  in  Sections  3.3  through  3.5. 

f 

3.3.  Determination  of  H,  E,  and  D 

!  Since  the  plane  is  subjected  to  uniform  tension,  p, 

•  1 

[.arallel  to  the  axis  of  x  at  infinity  (see  Fig.  1)  ,  the  regu¬ 
larity  condition  (2.23)  is  then  sacisfied  at  infinity.  This 
can  also  be  expressed  in  terms  of  stress  function  as 

;  ' 

;  See  Equations  (2.20)  through  (2.22). 


2  2 

U  =  2^—  =  2  c ^  sinh2£  sin2n  =  (e2^+  e”^  -  2)  (1  -  cos2n) 


Since  at  £  both  e  ^  and  2  can  be  neglected  compared  with 

2£ 

e  ,  U  becomes  at  £  00 


q  =  e25  -  E|!  a25  cos2n 

16  16 


(3.22) 


From  (3.1)  and  (3.4)  through  (3.7),  U  becomes  at  £  -*•  °° 


U  =  He2^  +  Ee2^  cos2n 


(3.23) 


Comparing  (3.22)  with  (3.23),  w<s  get 


h  -  2£_ 
H  16 


(3.24) 


E  =  -  2£_ 

16 


(3.25) 


Notice  that  D  is  the  coefficient  of  U^.  Since  U4  is  inde¬ 
pendent  of  couple-stress  effects,  D  is  assumed  to  be  in  the 

"ic 

same  form  as  one  obtained  in  the  classical  case  [8],  namely, 


D  =  -2|~(cosh2£Q  -  1) 


(3.26) 


3.4.  Determination  of  the  Relation  Between  B,  F,  and  p  by 
Use  of  Cauchy-Riemann  Equation's,  (2.13) 

Using  the  expressions  in  (3.9)  through  (3.15)  together 

with  (3.6)  and  (3.7),  we  have 


See  Equation  (6.2495)  on  p.  544. 
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(3.27) 


v2u  =  hv2u1  +  bv2u2 


T  -  £2V24'  =  F(4’1  - 


£2V24'1)  =  F^ 


(3.28) 


By  direct  substitution  of  (3.27)  and  (3.28)  into  (2.13),  the 
Cauchy-Riemann  equations  read 

F  3^  =  -2(1  -  v)Jl2[H~(V2ui)  +  b|^(V2U2)]  (3.29) 

3  4* 

F  — p  =2(1-  v)a2[H^-(V2U1)  +  b|^-(V2  U2)]  (3.30) 

Using  the  results  obtained  in  (3.9),  (3.10),  (3.14),  and 
(3.24),  either  (3.29)  or  (3.30)  leads  to  the  same  expression 

—£=  ( 1  -  v )  £2  (B 1  j  -  p)  (3.31) 

c  c 


or 


B 


Fc 


16(1  -  v)4‘ 


(3.32) 


3.5.  Determination  of  C,  F,  and  M 

The  traction  conditions  on  the  interior  boundary  will 
now  be  used.  On  the  boundary,  £  =  and  g  varies  within  the 
limit  (0,2tt)  . 

a.  Couple-stress  On  the  boundary,  expression 
(3.20)  becomes 
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•V<5o'nl  = 


cA 


✓7 

j^/2'  ~2F  smh2^o  sm2n  + 


Ml(Fek-)  ,)r 
2 


•  [ (j  +  cos^2^q) (se2)  -  (2  cosh2?o)  (se2>*  (3.33) 

•  cos2ri  +  ^-(se2)  cos4n]  r 

The  boundary  condition  (2.22)  is  now  applied  by  multiplying 
(3.33)  on  both  sides  with  se2,  and  by  integrating  the  result¬ 
ing  expression  with  respect  to  n  through  the  range  (0,2tt). 

The  following  expression  results  with  the  use  of  equations 
developed  in  Appendix  II  (AII.l)  at  n  =  2 ,  (All. 4)  and 
(All. 5)  at  n  =  0,  p  =  1,2. 

M  =  2TF  (3.34) 


v;here 


1’  =  B~2^  sinh2£  /  (Fek„ )  '  ]  ,  ( (i  +  cosh22£  ) 

4  O  /  s  4  O 

O 


+  (2 


cosh2£Q) 


I 

r=0 


B 


(2) 

2r+2 


B(2)  + 

B2r+4 


I 

r=0 


B 


(2)  , 

2r+2 


(3.35) 


i(2)  | 
2r+6! 


2 

)  ]} 


b.  a£ri  on  the  boundary.  After  (3.18)  is  evaluated  on 
the  boundary,  further  reduction  is  introduced  by  the  use  of 
(3.24)  ,  (3.25) ,  (3.26)  ,  and  (3.32)  . 
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pjfU' 


+  c°sh^25o)  -  £2(cosh2£Q)  +  t^Jsii^n 

^"Co 


+  [-?! (cosh2SQ)  +  C2(y  +  cosh22CQ)  +  ?4]  sin4n 


1  O 

+  f~3 —  ^2  cosh2?o]  sin6n  +  [_£]sin8n 


~  [55(5-  +  cosh22€Q)  cosh2£;o]  (se2)" 


+  [?5(3  cosh22CG  +  j)  ]  (se2)"  cos2n 


(3.36) 


-  t^5(4  cosh2CQ)]  (se2)"  cos4n  -r  [-1]  (se2)"  cos6n 

-  (Sg(*r  +  cosh22C0)j  (se2)  +  [5g(2  cosh2£Q)]* 


(se2)  cos2n 


-  [-J-]  (se2)  cos4n  +  [£7(  j  +  cosh22£Q  )]  (se2)  '  sin2n 


^  7 

-  U?  (cosh2CQ)  ]  (se2)  '  sin4n  +  [-4](se2)'  sin6n 


where  K  is  a  constant  and 


=  - 


8(1  -  v)r 


-4E 

F  cosh2E  -  C(e  0  +  3) 


2  2  4£, 

-  ^~(cosh2f;o  -  1)  -  E^— ( e  0  +  3) 
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t?  =  2Fr(Fek2)c=c 

The  boundary  condition  (2.21)  is  used  by  multiplying  (3.36) 
on  both  sides  with  se2  and  by  integrating  with  respect  to  n 
from  0  to  2tt.  Further  simplification  is  achieved  by  using 
(AII.l)  with  n  =  2, 4, 6, 8;  (All. 4)  and  (All. 5)  with  n  =  0, 
p  =  1/2;  (All. 6)  with  n  =  0,  p  =  1,2,3;  and  (All. 7)  with 
n  =  0,  p  =  1,2,3.  All  algebraic  manipulations  are  routine 
but  cumbersome,  thus  omitted.  The  resulting  expression  is 

bC  +  gF  =  m  (3.37) 

where 

b  =  j-[(e  4?°  +  3)  (j  +  cosh22SQ)  +  e  ?°  cosh2CQ]  B^2) 

_4  £  _2£ 

[  (  e  °  +  3)  (cosh2£0)  +  e  °( £  +  cosh22?o)]  B^2) 
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**!■-*  r-i—,  *M! 
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mswmmmmrn 


2  4£ 

m  =  J|_{ [ (e  °  +  4  cosh25( 


3)  ( +  cosh225  > 


2 5rt  ,,, 

+  e  cosh.,5  J  B ' 
o  2 


+  l  (-e  °  +  4  cosh25  **  D  cosh25 

o  o 

+  e  ^°(  J  +  cosh225Q)lB^2) 


(3.40) 


r 1  .  4S0 

+  lj(e 


+  4  cosh2E  -  1)  +  e 
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°  cesh25  ]Bi2) 

O  b 


+ 


B 


c.  on  the  boundary.  Equation  (3.16)  is  evaluated 

first  or.  the  boundary.  The  following  expression  results  with 
the  substitution  of  Equations  <3.24;;,  (3.25),  (3.26),  and 
(3.32) . 


[ K'O, J 

”  5=5. 


y 

l"T  +  Y2  cosh25Q  +  +  coi;h2250)  -  Y4J 

2  3 

+  t -y,  cosh25Q  -  y2 (cosh  25Q  +  j) 

-  Y->(2  cosh25.J  +  yJ  cos2n 

J  o  o 

1  2 

+  t  +  cosh  250)  +  Y2(cosh250) 


+ 


Y 

+  y5]  cos4n 
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+  1-y^  cosh2?o  -  -y-]  cos6n  +  [-y^cos8n 


(3-41) 


Iy*,^  +  cosh22£  )  cosh2€ 
7  2  O  < 


-  Yg(^  +  cosh?2Cd)l  (se2) ' 


+  (y^(3  cosh225Q  + 


-  y g  ( 2  cosh2?o)l  •  (se2) '  cos2n 


-  f  \  coshlK)  — (se2)  '  cos4n 


y7 

+  [  -^-1  (se2) '  cos6n 
+  l  Yg  +  cosh22£0)]  •  (se2)  sin2n 


-  [  Yg  cosh2?oI  (se2)  sin4n  +  (-j-Jfee^)  sin6n  =  0 


where  K'  is  a  constant  and 


?,  „  2f'c*  .  Fc2 


rl  =  ff(2  "  e  ^  + 


+  Ce 


16(1  -  v>r 


1  2  2  4E 

y  =  - — - -  cosh2E  +  f|-(8  cosh2E  -  e  -  3) 

2  4(1  -  v)£2  °  ^ 


+  (e  +  3)  C 
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r 


JS—nWi.  > '•■c-y.'-y  w&zrr**’*'  mwiw1* 


~  2  -4£  2  4£  -4£ 

V,  -  - ^ - j,2  +  e  °)  ♦  |f-<4  +  e  °  4  e  ° 


16(1  -  v) £7 


2 p  -2^ 

-  3e  °)  +  3Ce  °  +  D  r.inh2£- 


y-  -  9F  sinh2£ 
4  < 


Yc  =  3F  sinh2£ 

5  o 


Y r  -  3F  sinh4? 

6  O 


Y?  =  2rF[(Fek2) ']?=5 


y8  -  2rF[(Fek2)']w 


Y9  =  2rF(Fek2)^_^  sinh2C( 


The  boundary  corlition  (2.20)  is  applied  with  the  multiplica¬ 
tion  of  (3.^1)  by  ce2,  then  by  integrating  the  resulting 
expression  with  respect  to  n  from  0  to  2ir.  By  using  (All. 2) 
and  (All. 3)  with  n  =  2,  4,  6,  8;  (All. 8)  with  n  =  0,  p  =  1,2,3; 
and  (All. 9)  with  n  =  0,  p  =  1,2,3,  the  final  expression  re¬ 
sults  in 


kC  +  sF  =  t 


(3.42) 


where 
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7 


k 


wqra mtmmmm  **"G* ■  ■  mtev- ' 


7  -2Crt  -4£  -2C 

<  — {-Ij  e  +  (2e  +6)  cosh2£Q  +  6e  ° 


cosh225  JA^2^ 

Q  O 


-25 

-  [7e  °  cosh25 


-45 

+  (e  0  +  3)  •  (cosh2250  +  |)Ja]2) 


(3.43) 


1 2e  0  +  e  0  cosh22£  +  (e  0  +  3)  *  cosh2£  ]A.(?) 


’oJ  4 


[e  0  cosh25  +  j(e  +  3)]A*2)  -  ^ — A*2)} 
0  4  6  4  8 


-  r  c 

0—1  ’ 


45, 


16(1  -  v) £ 


2<-Ij+8  cosh225Q  +  (2  +  e  °)  (1+2  cosh225Q} 


(18  sinh25Q/(  — ^iii^)]A^2) 

c 


[cosh25o  +  4(cosh25o)(|  +  cosh225Q) 


-45 

+  (2  +  e  °) (2  cosh25  ) 

o 


,3  si„h4{0)(iiil^jW»:nAW 


12  7  1  —45 

-  [(^  +  cosh  25  )  +  4  cosh  25  +  (1  +  i-  e  °) 

^  o  o  z, 


+  (3  smh25o)(  16(1  -v)l  )]AU) 

c 
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[2  cosh2C0JA^2)  -  (j)a|2)} 


-1  9 

2r(  [  (=•  *  cosh  2f  )  (sxnh2£  )  (Fek_) 


(|  +  cosh22CQ)  (Fek2) •  cosh2£  ] 


HI 

I.  <2r  +  2)?.^!,  B<2> 


L  r=0 


2r+2  2r+2 ' 


+  [2(cosh2CQ) (sinh2^o) (Fek-) 

S  \>r 


-  (3  cosh22^o  +  jHFek2)^=^  ]  •  [a^2)B^2) 


*7  1  <2r  +  2)a‘2»4  b(2|2] 


r=0 


+  [|  sinh2^&(Fek2)?=?  -  |(cosh2£0)  (Fek2U=^  ] 


•  rAf>B<2>  *  2A^2>Bj2)  -  2  f  (2r  +  2)A<2i,  B>"  J  ) 

r=0 


(2) 

2r+4  “2r+2J 


+  2r(Fek2)?=?  {-  i[A^2)B2(2)  +  2Aj2)B^2) 


+  3A^2)B^2>  +  i  f  (2r  +  2]A<2i0  B<?ij 


r=0 


2r+8  2r+2J 
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+  III  +  cosh22£  |  r  ^  A{2)  B(2)  -  y  A(2)B(2} 
2'4  S'  -  A2r+2  2r+2  A2r  2r+2 

,r=0  r=0 


-  —  (cosh25  )  Y  A^)  B^2^  -  s  A^2^B^2^  ' 
2(cosn^^o;^^  A2r+6  B2r^-2  r^Q  A2r  2r+6J 


if  *  a(2)  (2)  _  r  a(2)b(2)  l 

2>  A2r+8  B2r+2  A2r  B2r+8  / 

Lr=0  r=0  J/ 


(3.44) 


2  -  25  * 

?§-U  i  (e  °-2)  +  2(3  +  e 


16  i  2 


8  cosh25  )  cosh25 
o  o 


+  (e25°  +  2e  2?0  -  2e  4?0  -  4) (1  +  2  cosh225Q) ] A^2) 


2?  -25  -45 

+  t  (-3e  °  -  4e  °  +  4e  °  +  10)  cosh25 


-  (3  +  e  °  -  8  cosh25Q)  (cosh225Q  +  |)]A^2) 


-45  4  5 

-  t(2  cosh25  -  e  °  -  3)  +  (e  °  +  3)cosh25 


25  o  ,n\  25  ,  ,  4E 

+  (e  °  -  10)  cosh225  ]a;2)  -  [e  cosh2£  +  j  +  ~  e 

O  4  0  4  4 


4  cosh2S  ]A(2)  -  i(e  ’°  -  2)A^2) 

Ob4  o 


(3.45) 


Equations  (3.37)  and  (3.42)  form  two  linear  nonhomogene- 
ous  equations  with  two  unknowns,  C  and  F.  Xhey  can  readily  be 


solved  by  using  Cramer's  rule 
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bs  -  gic 


(3.46) 


_  bt  -  ni< 
bs  -  q< 


(3.47) 


provided  (bs  -  g<)  =  0. 

The  stresses  and  couple-stresses  of  this  problem  are  the 
general  expressions  in  (3.16)  through  (3.21),  together  with 
the  seven  determined  coefficients  in  (3.24)  through  (3.26), 
(3.32),  (3.34),  (3.46),  and  (3.47). 


4.0  SPECIAL  CASES  OF  THE  S 


10 


4.1. 


Infinite  Plane  with  a  C 
Tension  in  the  Presence 


o 

e 


Rewrite  (3.6)  and  (3.7) 


U  =  -^  (c2Ux)  +  (B) (U, 

c  1 


+  (D)  (U4)  +  fc" 

c 


<l>  =  (F) 


0 


+  (M) 


IT 

?  r\  \ 


Solution  of  this  specific  problem  can  be  obtained  by  use  of  the 
ximiting  process  described  in  Appendix  III- (a).  Equations  (4.1) 
and  (4.2)  are  simplified  and  reduced,  by  use  of  the  forms  in 
(AIII.6)  and  (AIII.7)  and  the  basic  theorem  on  limits 


U 


=  j  r^(l  -  cos  <J>)  -  £n  r 


(4.3) 


2T?V 


AT 


a  (1  -  A) 

2p 


+  1 


cos  2 <t> 


pafc 


1  +  A 


4(1  -  v)a£ 


sin  2<p 


(4.4) 


which  are  identical  to  the  forms  obtained  by  [2].* 


4.2.  Infinite  Plane  with  a  Crack  under  Simple  Uniform  Tension 
(Fig.  6) 

Solution  to  this  problem  can  be  obtained  from  the  solu¬ 
tion  of  the  elliptical  hole  problem  by  taking  ■+■  0,  as 
described  in  Appendix  III-(b).  The  results  in  (AIII.19)  show 
that  in  this  crack  problem  the  couple-stress  effect  vanishes. 
Furthermore,  (AIII.19)  shows  that  the  solutions  obtained  here 
are  identical  to  the  corresponding  classical  solutions  which 
can  also  be  obtained  by  proceeding  by  the  sam®  limiting  process 
in  [8) .  Since  couple-stresses  are  related  to  curvature  and  no 
curvature  is  induced  in  this  crack  problem,  it  is  clear  that 


There  is  a  misprint  for  the  expression  of  E  in  Equation 
(38)  in  [2].  The  correct  form  is  E  =  (-pa£F)/(l  +  F)K^(a/£). 
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couple-stress  solutions  should  be  the  same  as  classical 
solutions. 


£q+o,  a*c 
b+o,  e+1 


Figure  6. 

4.3.  Classical  Solution  with  Free  Stresses  on  5  =  5Q  (Fig.  7) 
As  i  ■*  0,  couple-stress  solutions  of  the  problem  in 
Fig.  1  obtained  in  Section  3.0  reduce  to  classical  solutions 
of  the  problem  in  Fig.  7.  The  results  in  Equation  (AIV.9)  are 
identical  to  the  forms  obtained  in  [8].* 


pp. 


T 


See  Equations 
543-544  of  [8] . 


(6.2481) , 


(6.2494) 


and  (6.2495)  on 
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5.0  DISCUSSION  AND  CONCLUSION 


The  solutions  to  this  problem  are  given  either  as  stress 
and  couple -stress  functions  in  (3.6)  and  (3.7)  or  as  stresses 
and  couple-stresses  in  (3.16)  through  (3.21),  together  with 
the  seven  determined  coefficients  in  (3.24),  (3.25),  (3.26), 
(3.32),  (3.34),  (3.46),  and  (3.47).  All  series  in  (3.32), 
(3.24),  (3.46),  and  (3.47)  are  shown  to  be  convergent  in 
Appendix  IV.  Special  cases  of  the  solutions  are  also  obtained 
by  the  proper  limiting  process  as  discussed  in  Section  4.0. 

The  three  limiting  cases  are:  (4.3)  and  (4.4)  for  the  degen¬ 
erate  circle  in  Fig.  5,  (AIII.19)  for  the  degenerate  crack  in 
Fig.  6,  and  (AIV.9)  for  the  classical  solution  in  Fig.  7. 

The  result  obtained  is  a  subclass  of  solutions  to  the 
general  self-equilibrated  boundary-value  problems.  The  selec¬ 
tion  of  the  Mathieus'  functions  and  the  form  of  weighting 
functions  in  the  boundary  conditions  (2.20),  (2.21),  and 
(2.22)  match  a  particular  class  of  boundary-values,  which 
gratifyingly  does  reduce  upon  limiting  processes  to  various 
classical  solutions  (Figs.  6  and  7,  and  Fig.  5  for  1=0). 

Here  we  summarize  some  important  results  about  the 
three  limiting  cases.  The  solutions  for  the  degenerate  circle 
in  (4.3)  and  (4.4)  are  shown  to  be  identical  to  the  results  in 
[2].  The  solutions  for  the  degenerate  crack  m  (AIII.19)  are 
seen  to  be  the  same  as  those  obtained  in  [8],  by  proceeding 
according  to  the  same  limiting  process.  Since  couple-stresses 
are  related  to  curvature  and  no  curvature  is  induced  in  this 
crack  problem,  it  is  clear  that  couple-stress  solutions  in 
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(AIII.19)  should  be  the  same  as  the  classical  solutions  in  [8]. 


The  classical  solutions  for  the  general  elliptical 
lem  (Fig.  7)  in  (AIV.9)  are  obtained  by  taking  £  -* 
and  (3.7)  and  they  are  identical  to  the  results  in 


hole  prob- 
0  from  (3.6) 
18]. 


1 


APPENDIX  I 

Solution  of  a  Two-dimensional  Wave  Equation  in  Elliptical 
Coordinates 

The  two-dimensional  wave  equation 
2  2 

7  -  £  v  7=  0 


can  be  expressed  in  elliptical  coordinates  £  and  n  by  using 
(2.12)  as 


+  1*  _  2q  (cosh2£  -  cos2n) 7  =  0  (AI.l) 

3r  3n 

Using  the  method  of  separation  of  variables  and  the  form  of 
solution  of  (AI.l)  as  7(5, n)  =  7(5)<Mn)»  (AI.l)  is  reduced  to 
the  following  two  ordinary  differential  equations  known  as 
Mathieu's  equations. 


2 

d  7 


(d  +  2q  cosh25)7  =  0 


(d  +  2q  cos2n)  =  0 


(AI.2) 


(AI.3) 


Here  d  is  the  separation  constant. 

From  the  solutions  of  (AI.2),*  we  choose  the  following 

one 


See  Section  8.30  on  p.  165  and  Section  13.30  [Equation 
(5)  on  p.  248]  in  Reference  [9]. 
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Fek2(C/-q)  = 


(P2) ’  “  (2) 

—ToT  L  KZ  I  (v.JK  (v9) 
irA^2*  r=0  2r  r  1  r  2 


(AI.4) 


From  the  solutions  of  (AI.3),  we  select,  for  the  current 
problem , 


se9(n,-q)  =  Y.  (” 1) r® wlo  sin  (2r  +  2)n 
z  r=0 


(AI.5) 


However,  another  Mathieu  function  of  (AI.3) 


ce9  (  n,~q)  = 


00 

Y_  f-l)rA92^  cos2rn 


(AI.6) 


will  be  used  in  the  text  for  evaluation  of  integrals. 


Derivatives  of  se2(n,-q),  Fek2(£,~q) 

Since  the  series  expressed  in  (AI.4)  and  (AI.5)  are 
uniformly  convergent**  for  real  £  and  n>  respectively,  they 
may  be  differentiated  term  by  term. 


^•[Fek2(£,-q)  ]  = 


(p2) 'k  « 


—ITT  I  le^1  (v.  ) — 5— 

ttA(2)  r=0  2r  r  1  3v 
o 


8Kr(v2) 


_r  1) 

“  e  W] 


See  Section  2.18  on  p.  21  in  Reference  [9]. 


(AI.7) 


See  Sections  3.21  and  3.22  (on  pp.  37-38)  and  Section 
13.60  on  p.  257  in  Reference  [9). 
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(AI.8) 


g^-lse2(n,-q)]  =  £  (-l)r(2r  +  2)B^2{2  cos(2r  +  2)n 


^[se^m-q))  =  I_  (-Dr+1(2r  +  2)2B^|2  sin(2r+2)n  (AI.9) 
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'  *3"' 


■~r»Ka f*w 
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lse2n+2(n,-q) J  cos2  pndn  = 


(P-3) 

y  E(2n+2)p(2n+2)  _ 

r=0  2r+2  2p-2r-2 


_  y  R(2n+2)  R(2n+2) 
yzQ  °2r+2  °2p+2r+2 


(p  is  odd,  >1)  >  (All. 5) 


¥ti2a+2y 


(P-4) 

2  R(2n+2)  _(2n+2)  . 
£  B2r+2  B2p-2r-2 
r=0 


V  R (2n+2)  (2n+2) 

rt0  2r+2  B2p+2r+2 


(p  is  even,  >  2) 


1  f2lf 

^  J  fse2n+2(T1'-q)1  fse2n+2(T1/-q)  J  sin2  pndn  “1 


/_nP  n  V  R (2n+2)  (2n+2) 

P  rt0  2r+2  2p+2r+2  (p  1  1} 


(All. 6) 


se2n+2(l''-q)se2n+2(,,'-'3)  = 


l  (2r+2)B<^2> 
r=0 


i  r2ir 

”  J  se2n+2(n,"q)se2n+2(n ,_q)  cos2Pndn  = 


(All. 7) 
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/ 


!<P-3> 


-  I 

r=0 


(2p  -  2r  -  3) 2 


B (2n+2) 
®2r+2 


„(2n+2) 

2p-2r-2 


(All.  7) 


+ 


00 


I 

r=0 


(2p  +  2r  +  l)2 


(2n+2) 

2r+2 


(2n+2) 

2p+2r+2 


(p  is  odd, >  1) 


< 

-  i  PB(2n+2) 

2  *p 


(p-4) 

2 


lo 


(2p  -  2r  -  3) 2  B^22> 


_ (2n+2) 
B2p-2r-2 


V 


I 


r=0 


(2p  +  2r  +  1)  B 


(2n+2) 

2r+2 


(2n+2) 

2p+2r+2 


(p  is  even,  >  2) 


2  r2iT 

?  JQ  Ce2n+2tn'~q)se2n+2*n'-q) 


sin2pndri  = 


\ 


T  k[2nt2)  3(2n+2)  -  V 


r=0 


2r+2p  2r+2 


(2n+2)  (2n+2) 

L.  A2r+2 
r=0  2r 2 


2r+2p 


)  (All. 8) 


(P  >  1) 


i  f2" 

*  -*0 


ce2n+2(n'-‘!)se2n+2(T|''<J)dn  = 


I 

r=0 


(2r  +  2)  a(2n+2)  (2n+2) 

1  2  2r+2  B2r+2 


2 

TT 


2ir 


ce2n+2  (n,_q)se2n+2  (n,-q)  cos2  pndn  = 


'2n+2 


J(P-I) 

V 

L. 

r=0 


r=0 


B2r+2 

+ 

(All. 9) 

(2n+2) 

(p  is  odd,  >  1) 

2r+2 
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II  M  8 
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APPENDIX  III 

Let  r  and  r'  represent  the  semi -major  axis  and  the  semi¬ 
minor  axis,  respectively,  for  any  confocal  ellipse.  By  using 
(2.8)  r,r'  can  be  expressed  in  terms  of  c  and  K 

r  =  x  =  ±c  cosh?,  at  n  =  0,ir,  y  =  0 

(AIII.l) 

r'  =  y  =  ±c  sinh?,  at  n  =  ^-rr,  x  =  0 

By  definition 

c  =  re  (AI1I.2) 

From  (AIII.l)  and  (AIII.2),  we  have 

cosh?  =  e"1  (AIII.3) 

Some  Useful  Limiting  Forms  When  the  Ellipse  Tends  to  a  Circle 
with  e  -*■  0 

For  fixed  r,  as  e  +  0,  the  confocal  ellipse  of  the 
semi-major  axis  r,  tends  to  a  circle  with  r  as  radius,  and  the 
confocal  hyperbolas  become  radii  of  the  circle,  with  n  =  <}> 

(see  Figure  3(b)).  By  use  of  (AIJI.2)  and  (AIII.3),  it  is 
clear  that  e  ■*  0  com:,  sponds  to  c  ■*  0  and  ?-*•«>.  Since 

=  k  =  i  £  .  <  •'  j  or  k  0  as  c  0.  By  means  of  these 

/  X/ 

•ff 

equivalent  limiting  processes,  we  obtain  the  following 

lim(i  ce^)±n  =  rtn  (n  =  1,2,3)  (AIII.4) 

e-0 

"jf 

See  Appendix  I  on  p.  367-370  ir  Reference  [9] . 
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(AIII.5) 


lim  7FTr  Fek2l 

e-*0 


lim  7~~c~r  (Fek„)  1 

e+ol^  2 


r 


3r 


K, 


(AIII.6; 


lim(se-)  =  sin2<(>  (AIII.7) 

e-*-0 


lim(ce-)  =  cos2<j>  (AIH.8) 

e-0  z 


For  the  specific  ellipse  of  the  semi-major  axis  a  with  £  =  £  , 
(AIII.7)  and  (AIII.8)  still  hold  while  (AIII.4)  through 
(AIII.6)  take  the  following  forms. 


lim 

e-*0 


’  1 

±n 

[2  ce 

±n 


(n  =  1,2,3, ... )  (AIII.9) 


lim 

e-*0 


TT 


(Fek2) 


(AIII.10) 


lim 

e-*-0 


(Fek2) 


(AIII.ll) 


Some  Useful  Limiting  Forms  When  the  Ellipse  Tends  to  a  Crack 
with  e  +  1  or  4  0. 

The  ellipse  with  the  semi-majo.r  axis,  a  =  c  cosh  £  ,  and 
the  semi-minor  a,. is,  b  =  c  sinh£  ,  tends  to  the  inter-focal 
line  of  length  2c,  when  e  -*■  1,  while  a  ->•  c,  b  -*■  0  (Figure  3(b)). 
By  use  of  (All I. 3)  with  C  =  CQ,  e  ^  1  is  seen  to  be  equivalent 
to  CQ  4  0.  The  following  useful  limiting  quantities  are 
obtained  by  direct  computation. 
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lim  e 


(n  =  1.2,3,...) 


(AIII.12) 


lim[smh2£  ]  =  0 


lim[cosh2£  ]  =  1 

V° 


(AIII.13) 


lirai(Fek_)  i  =  Fek_(0,-q) 

£  -0  ^ 


(AIII.14) 


lira  (Fek„)  r  ,  j  =  Fek  '0,-q)  = 

C  -*-0  I  ^  V  ^ 


P2ce2 


(AITI.15) 


Limiting  Values  of  Solutions  as  e  0 

By  means  of  the  forms  in  (AIII.4)  through  (AIII.ll),  the 
limiting  values  of  the  modified  forms  of  (i  =  1,2,.,., 5)  in 
(3.1)  and  the  modified  forms  of  4,^,4,2  in  (3.4)  and  (3.5)  are 
obtained  as  follows 


lim[c2U1]  ?  4r2 
e-*-0 


lim[U-]  =  cos2<}) 
e+0 


lim  --S 
e+Olc  - 


cos<J> 


lim[U. ] 
e+0 


=  In  r 


(AIII .16) 


lim[c2Uc]  =  4r2  cos2<)> 
e-0  5 
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e+0  c 


~j  sin2<j> 


i“[i5^r  v2j  -  k2(i)  sin2» 


By  use  of  the  forms  in  (AIII.4)  through  (AIII.ll),  the  limit- 

2  * 

ing  values  of  D  in  (3.26),  B  in  (3.32),  c  C  in  (3.46),  F  in 
**  +  2 

(3.47),  '  T/c  in  (3.35),  and  M  in  (3.34)  are  obtained  as 

follows. 


lim  [D]  =  - 
e+0 


lim[d  C]  = 
e+0 


>a  (1  -  A) 
(1  +  A) 


lim[F]  = 
e-*-0 


4(1  -  v)  a2£2f 

— rr+ at 


lim[B] 

e->-0 


2 (1  +  A) 


(AIII.17) 


- 


(p2}  16a&  (1  -  v)^  - 


c2d  is  obtained  by  multiplying  (3.46)  by  c2. 

*  *  ~ 

F  is  obtained  by  multiplying  (3.47)  by  cz/2. 

^b,  g,  M,  n,  s,  and  t  in  (3.46)  are  expressed  in  (3.38), 
[3.39),  (3.40),  (3.43),  (3.44),  and  (3.45),  respectively. 

t  2  9 

T/c  is  obtained  by  dividing  (3.35)  by  c^. 
is  obtained  by  multiplying  (3.34)  by  (p2) '/1T* 
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e+0 


(1  +  v)^  J 


where 


8(1  -  v) 

K  - 

a  ^  a  o  i 
I  +  2  I  ~T 
Ki  z 


2  2 

Since,  from  (3.24)  and  (3.25),  we  have  H/c  =  -  E/c  =  p/16 
which  is  a  constant,  the  following  limiting  processes  are 
obvious . 


■  i  •  H 

lim  — ~- 
e+0  c 


e+0  c 


-  JL 

16 


(AIII.18) 


Limiting  Values  of  Solutions  as  0 

By  means  of  the  forms  in  (AIII.12)  through  (AIII.15), 
the  limiting  values  of  H  in  (3.24),  B  in  (3.32),  C  in  (3.46), 
D  in  (3.26),  E  in  (3.25),  F  in  (3.47),  and  M  in  (3.34)  can  be 
obtained  as  follows. 


lim  H  =  lim  b  = 


V° 


V° 


lim  C  =  lim  E  = 


v° 


v° 


-  E£i 
16 


(AIT 1.13) 
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APPENDIX  IV 


Some  Useful  Limiting  Forms  as  &  -»  P 
1  /2  1  c  ' 

Since  q  /  =  -j;-*-  =  k,  for  fixed  c,  0  implies  that 

q  4  ®  or  k  +  The  following  limiting  quam  .ties  are 
obtained  by  using  [9].* 


A<2) 

lim  — T^r-  =  2  (-l)r  (r  =  1,2,3,...)  (AIV.1) 

A<2) 


b(2) 

lim  -%±2  =  (-l)r  (r  +  1)  (r  =  1,2,3,...)  (AIV.2) 

*■*•0  B^  ' 


lim 


=  0 


(r  =  0,1,2, 3, .. .)  (AIV.3) 


lim 

£-*0 


0 


(r  =  1,2,3, . . . )  (AIV.4) 


lim 

£->0 


lim 

£->-0 


1 


(AIV.  5) 


For  large  arguments, 


*  * 

z ,  we  have 


I  (z) 
n 


I'  (z) 
n 


(AIV. 6) 


-z 

K  (z)  — - - ~ K  '  ( Z )  (AIV.  7) 

n  /2x7^ 


s - 

See  Section  3.34  on  p.  47. 

* *See  Equations  (9.7.1)  to  (9.7.4)  on  pp.  377-378  in  [9], 
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Substituting  (AIII.6)  and  (AIII.7)  into  (AI.4)  and  (AI.7)  for 
the  large  value  of  q  (small  value  of  £) ,  and  taking  the  limit 
of  the  resulting  form  with  £  -*■  0,  we  obtain 


lim 

£->0 


[Fek2] 

[Fek2] 


«o 


sinh£ 

o 

2 


(AIV.8) 


All  the  series  forms  in  expression  g  in  (3.39)  and  s  in  (3.44) 
can  be  shown  to  be  absolutely  and  uniformly  convergent.*  As 
such,  by  use  of  [11],**  (AIV.3)  and  (AIV.4) ,  all  these  series 
tend  to  be  zero  as  £  0.  Based  upon  this  argument  and 

(AIV.8),  it  is  seen  that  in  g/B2  and  s/A^  ,  terms  contain¬ 
ing  £  as  denominators  dominate  for  the  small  value  of  £. 

(2)  12) 

Hence,  other  terms  in  g/B2  and  s/A^  can  be  neglected  as 
compared  with  terms  havinq  £  as  denominators,  for  £  ■*  0. 

Limiting  Values  of  Solutions  Like  £  ->  0 

By  use  of  (AIV.l)  through  (AIV.8),  the  limiting  value 
of  H  in  (3.24),  B  in  (3.32),  C  in  (3.46),  D  in  (3.26),  E  in 
(3.25),  F  in  (3.47),  and  M  in  (3.34)  can  be  obtained  as 
follows . 


lim  [H] 
£->0 


1 

16 


pc 


2 


1  2 

lim  [B]  =  yv-  Pc 
£-*-0  ib 

5 

See  Section  3.22  on  p.  38  in  [9], 

ft  ft 

See  Theorem  2  on  pp.  339-3  iO. 
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